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INTRODUCTION 


In  recent  years  military  equipment  designers  have  become 
more  and  more  interested  in  the  behavior  of  dust  particles  in  an  air 
stream.  The  military  environment  is  often  a  dusty  one,  and  dust  is 
well  known  as  a  destroyer  of  combustion  engines.  Several  general 
types  of  air  cleaners  have  been  used  to  clean  the  air  for  use  in 
powered  ground  equipment.  These  cleaners  have  usually  been  based  on 
inertial,  electrostatic,  or  impingement  principles,  and  their  design 
has  required  compromises,  particularly  in  the  areas  of  filtration 
efficiency,  pressure  drop,  size,  weight,  complexity,  and  service  life. 
The  advent  of  gas  turbine  engines  in  military  equipment,  increased 
use  of  aircraft  in  the  dusty  region  near  the  ground,  and  the  general 
push  toward  increased  efficiency  with  decreased  size  and  weight,  has 
brought  air  cleaner  designers  to  the  limits  of  available  basic 
knowledge.  One  of  the  most  lucrative  areas  of  ignorance  tempting  the 
researcher  in  fluid  dynamics  is  the  behavior  of  dust  particles  in  the 
boundary  layer.  Fluid  dynami cists  have  vigorously  attacked  the 
question  "what  happens  to  a  fluid  in  the  vicinity  of  a  wall,  where 
viscous  shearing  forces  are  relatively  large?"  Prandtl's  boundary 
layer  hypothesis,  and  later  work  by  such  able  researchers  as 
von  Karman,  have  permitted  the  aerodynamicist  to  work  in  an  atmosphere 
of  relative  certitude.  The  laminar  portions  of  the  boundary  layer  are 
susceptible  to  mathematical  solution,  although  the  procedures  are 
somewhat  sophisticated.  Turbulent  boundary  layers  are  less  amenable 
to  analysis,  although  statistical  techniques  offer  some  help  in  their 
description.  When  one  adds  solid  dust  particles  to  these  already- 
.  complex  flow  phenomena,  it  is  obvious  that  theoretical  analyses  wili 
be  difficult.  Of  equal  difficulty  are  the  techniques  needed  to 
measure  the  motion  of  the  small  dust  particles.  Optical  methods  used 
for  studying  larger  particle  movement  are  totally  inadequate  to 
observe  particles  of  the  order  of  one  micron  in  diameter. 
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The  vortex  separators  being  studied  by  Dr.  Hans  J.P. 
von  Ohain  and  his  associates  at  the  USAF  Aerospace  Research 
Laboratories,  the  paral lei -plate  dust  separators  studied  by  Dr.  Franz 
J.  Jaklitsch  at  ATAC,  and  the  improvements  on  impingement  filters  by 
Mr.  Richard  W.  Child  of  the  ATAC  Propulsion  Systems  Laboratory  and 
his  associates,  all  involve  motion  of  dust  and  air  in  the  boundary 
layer.  These  efforts  have  stimulated  the  present  study,  which  is 
pointed  specifically  at  the  geometrical  I y -simp test  device,  the 
Jaklitsch  paral lei -plate  filter. 

A  suitable  theoretical  development  has  been  achieved.  It 
is  based  on  the  earlier  work  of  the  Prandtl  and  von  Karman  schools, 
applying  additional  statistical  correlation  techniques.  For 
convenience,  the  development  assumes  spherical  particles  of  fixed 
density  compared  to  the  fluid  density*  There  does  not  anpear  to  be 
any  insurmountable  barrier  to  handling  variations  from  these 
assumptions.  It  is  also  assumed  that  particle  density  is  small 
enough  that  the  particles  behave  as  individuals,  with  negligible 
influence  from  other  particles.  This  assumption  is  basic  to  the 
theoretical  development  and  also  to  the  experimental  apparatus. 
Variation  from  this  assumption  appears  to  be  difficult  indeed,  and  is 
left  to  other  researchers.  Fortunately  there  appear  to  be  many  areas 
of  practical  interest  within  existing  assumptions.  The  instrumenta¬ 
tion  uses  a  carefully-collimated  light  beam  of  high  intensity, 
interrogated  at  right  angles  by  two  careful ly-col limated  optical 
seekers  of  known  separation  distance.  These  seekers  use  electronic 
light-amplification  techniques  to  provide  particle  count  and  particle 
transit-time  information  from  even  the  smallest  particles.  These 
data  are  collected  on  a  digital  tape  recorder.  The  raw  data  tapes 
are  subsequently  processed  in  a  general-purpose  digital  computer  to 
determine  the  desired  concentration  and  velocity  coefficients. 

It  has  become  apparent  that  the  instrumentation  system,  the 
Dobbs  Optical  Anemometer,  is  a  powerful  new  tool  for  research  in 
fluid  dynamics.  For  instance,  it  can,  by  using  very  small  particles 
as  tracers,  give  extremely  accurate  measurements  of  both  mainstream 
andboundary-layer  flow.  It  can,  in  an  expanded  form,  resolve. all 
velocity  components  of  interest,  including  turbulence  as  well  as 
average  velocities.  These  can  be  measured  with  no  appreciable 
instrumentation  influence  on  the  flow. 


t)  ,l"  ri**  •*  **'  •» 

0K  f.  MIHmH  (US  MT  N.  !•» 


198 


DOBBS  and  WISS 


ANALYTICAL  APPROACH 

A  1/4- Inch  wide  by  6-inch  high  slot,  32  inches  long,  set  in 
the  center  of  a  6-inch  square  duct  was  chosen  for  the  reasons  that: 

(1)  It  was  sufficiently  large  to  work  in;  (2)  It  was  long  enough  to 
Insure  fully  developed  turbulent  flow  profiles  in  the  slot;  (3)  The 
height-to-width  ratio, 24,  was  large  enough  to  permit  assumption  of 
two-dimensional  flow  at  the  center;  (4)  The  resulting  test  section 
was  small  enough  to 'be  built  without  undue  problems  (In  practice, 
this  has  proved  only  marginally  true*);  (5)  The  Reynolds  number  of  the 
flow,  based  on  the  hydraulic  diameter  of  the  slot,  would  be  high 
enough  for  turbulence. 

The  range  of  flow  velocities  chosen.  20  -  200  fps,  was  based 
on  the  results  of  the  earlier  work  at  ATAC  (2),  Lower  velocity  flows 
did  not  cause  deposition  of  particles,  and  higher  velocity  flows 
apparently  tended  to  strip  the  deposited  material  back  off  the  wails. 
The  resulting  range  of  variation  of  the  flow  Reynolds  number  is 
2500  to  25000. 

For  this  velocity  range  in  air  under  approximately  standard 
conditions  it  is  reasonable  to  assume  incompressible  flow.  Since  a 
filter  based  on  this  principle  would  normally  be  a  steady  flow  device, 
flow  was  assumed  time-steady  in  the  mean.  Finally,  a  particle  size 
range  of  0.1  to  10  micron  was  chosen.  Larger  particles  than  this 
would  be  more  influenced  by  gravitational  forces  than  by  fluid  turbu¬ 
lence.  Smaller  particles  present  additional  experimental  problems. 

In  general,  no  problem  in  turbulent  fluid  mechanics  has  yet 
been  solved  without  the  use  of  empirical  data  at  some  point.  The 
descriptive  equations  inevitably  contain  more  unknowns  than  the 
number  of  distinct  relationships  among  them.  Such  problems  have 
usually  been  attacked  by  assuming  that,  despite  the  lack  of  an 
applicable  uniqueness  theorem,  the  Navier-Stokes  equations  hold, 
instantaneously,  for  turbulent  flow.  This  assumption  is  implicitly 
made  here.  In  other  work  involving  discrete  particles  a  microscopic 
viewpoint  has  often  been  taken,  considering  the  behavior  of  the 
individual  particle.  In  all  cases  statistical  concepts  must  be 
applied  to  account  for  the  random  behavior  of  the  flow  and  the 
particles  in  it.  Here  the  concentration  by  weight  of  the  particulate 
material  is  assumed  to  be  low  enough  (below  1%)  so  that  the  behavior 
of  the  fluid  will  not  be  significantly  affected  by  the  particles. 

At  the  same  time,  due  to  their  small  size,  the  numerical  density  of 
particles  in  the  flow  (on  the  order  of  l(P/cc)  i$  high  enough  to 
permit  a  continuum  approach  to  the  variation  in  particle  density  in 
the  flow. 


The  most  applicable  previous  work  on  the  deposition  of 
particles  out  of  a  turbulent  flow  onto  parallel  bounding  surfaces 
was  done  by  Friedlander  (*0  in  1954.  He  worked  with  flow  through 
tubes  and  based  his  analysis  on  a  paper  by  Lin,  Moulton,  and 
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Putnam(5)  in  1 953  which  was  concerned  with  transfer  of  material  out 
of  solution  on  to  the  walls.  The  analytical  approach  in  both  cases 
was  to  divide  the  flow  into  a  viscous  sublayer,  a  turbulent  core,  and 
a  buffer  layer  between  them.  Eddy  diffusivity  was  assumed  to  occur, 
decreasing  in  effect,  all  the  way  to  the  wall.  The  main  stream  con¬ 
centration  was  assumed  to  remain  constant  in  the  direction  of  flow. 
Although  this  last  assumption  was  applicable  in  these  papers,  it  is 
precisely  the  change  in  concentration  in  the  direction  of  flow  that 
is  of  concern  in  this  paper. 

Applying  the  continuum  approach  to  the  particle  concentra¬ 
tion  in  the  flow,  consider  an  elementary  volume  of  space.  For  the 
x  direction  the  net  convective  change  is: 

£cUAyAz  +  (CUMyAzAxJ  -  CUdyaz  =  (CU)  4  y^zdx  (1) 

The  net  diffusive  change  is: 

[#t£  &  nwl-jffay**  =  7x(0T&)*yizix  (2) 

Summing  for  all  three  coordinate  directions  *  j  simplifying,  for  an 
incompressible  fluid  the  conservation  of  suspension  equation  is: 

■S*T  tcui>  +  TZjfefjii)  -  '  7$  (3) 

Assuming  instantaneous  concentration  and  velocity  can  be  expressed  as: 


C  =  C  +  c  (U) 

U  =  U  +  u  (5) 


and  using  Reynolds  rules  for  averaging,  the  turbulent  flow  form  of 
equation  (3)  is: 


(6) 


Consider  Figure  2.  For  a  slot  flow  U  s  U2  =  0.  Since  the  flow  is 
two-dimensional  and  there  is  no  change  along  the  z-coordinate, 
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Since  the  flow  is  time-steady  in  the  mean 


The  conservation  of  suspension  equation  for  a  turbulent  slot  flow  is 

^  +  v) +  b  (c“>  *  b  - 0  t7) 


II  III 


provided jtf,  the  molecular  diffusion  coefficient  is  assumed  constant. 

Equation  (7)  can  be  further  reduced  by  considering  the 
relative  magnitude  of  the  terms.  Experimental  data  on  slot  flow'®' 
indicates  term  IV  of  equation  (7)  to  be  no  larger  than  about  one 
per-cent  of  term  I.  Experimental  data  on  particle  deposition  rates'*1' 
indicates  that  term  I  is  small,  but  that  term  II  must  approach  zero. 
The  high  value  of  the  Schmidt  number  in  this  case  (nominally 
2.5  x  10^  to  5  x  1C)6)  implies  a  very  flat  concentration  profile 
falling  off  very  sharply  close  to  the  wall.  Terms  III  and  V  then 
should  have  high  values,  essentially  equal  to  each  other,  at  the 
edge  of  the  viscous  sublayer.  From  the  foregoing  it  may  be  concluded 
that  terms  I,  III,  and  V  of  equation  (7)  should  be  retained.  In 
dimensionless  form,  this  gives  the  following  approximate  equation 
for  the  particle  concentration  in  a  turbulent  slot  flows 

V  *£  -  JL_  £x,  +  is  =  0  (8) 

Jjr.  dq 

To  solve  this,  expressions  for  the  velocity,  V,  and  the  correlation 
function  g  are  needed.  Pai  found  a  solution  in  polynomial  form  for 
V.  This  is  (7)i  . 

V  =  A  +  B^z  +  C^n  (9) 


where  A,  B,  C,  n  are  constants  whose  exact  form  and  derivation  are 
not  of  concern  here. 

The  correlation  coefficient,  g,  must  be  measured  experi¬ 
mentally.  This  is  one  of  the  aims  of  this  work.  The  resulting 
curve  will  be  described  by  a  polynomial  in  ^  •  With  this,  any 
particular  solution  for  ^  ,  the  dimensionless  mean  concentration 
(no  matter  how  peculiar)  will  enable  the  following  transformation  to 
reduce  equation  (8)  to  a  homogenous  equation: 

¥  ■  ¥p  *  $  do) 

The  resulting  equation  is: 

-  R*Sc  (A  ♦  B r£  +  Cr|r*)  ”  0  0  0 

Separation  of  variables  reduces  this  to  two  ordinary  differential 
equations,  one  with  variables  coefficients,  which  may  be  solved  by 
,  ordinary  methods.  Substituting  for 
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7=  ecjp  i  (7  ) 

(12) 

These  are^  where  K  is  a  constant, 

6*  +  K0  =  0 

(13) 

4"  +  K  R*  Sc  (A  +  B-y2  +  C^n)/i  =  0 

0*0 

The  critical  points  in  this  procedure  are  the  determination 
of  an  expression  for  g  and  the  construction  of  a  matching  particular 
solution  for  ^ *  Without  experimental  data,  further  meaningful 
discussion  is  difficult  at  this  point. 

EXPERIMENTAL  APPROACH 

Many  physical  quantities  must  be  measured  in  an  experiment 
of  this  sort.  This  discussion,  however,  is  limited  to  those  involving 
the  particle  concentration  and  flow  velocity.  The  quantities  which 
must  be  measured  in  order  to  solve  and  evaluate  equation  (8)  are: 

(1)  the  mean  particle  concentration  at  a  series  of  points  across 
the  channel;  (2)  the  average  product  of  the  particle  concentration 
variation  and  the  y-direction  velocity  variation  at  a  similar  series 
of  points* 


Search  of  the  literature  disclosed  no  existing  method  of 
measuring  these.  The  discrete  nature  of  the  particles  indicated  the 
quantities  could  not  be  measured  directly  but  would  have  to  be 
computed  from  counts  taken  over  finite  space  and  time  intervals.  The 
problem  then  resolved  itself  into  one  of  how  to  restrict  these  inter¬ 
vals  to  small  enough  values  so  that  measurement  of  the  microstructure 
of  the  flow  would  not  be  averaged  out.  It  was  also  necessary,  due  to 
the  random  nature  of  the  data  to  be  observed,  that  the  measurement 
produce  the  least  possible  disturbance  in  the  flow. 

After  considering  and  discarding  a  number  of  counting  probe 
designs  it  was  decided  to  attempt  the  design  of  an  instrument  based 
on  light-scattering.  This  would  entirely  circumvent  the  problem  of 
disturbing  the  flow  when  measuring  it.  Light-scattering  particle 
counters  draw  the  mixture  of  fluid  and  particles  through  a  known, 
illuminated  volume  (normally  several  cc),  at  a  known  rate,  and 
count  all  particles  in  the  stream.  The  optical  anemometer,  as  the 
proposed  instrument  became  known,  would  not  have  this  control  over 
the  fluid  stream.  As' shown  in  Figure  3#  the  sensitive  volume  would 
be  defined  by  the  intersection  of  a  light  beam  and  a  field  of  view 
in  the  fluid  stream  itself.  This  sensitive  volume  had  to  be  kept 
as  small  as  possible  to  satisfy  the  basic  requirement  for  a  point 
measurement.  The  instrument  also  had  to  measure  the  "instantaneous" 
flow  velocity  at  the  point  where  the  particle  count  was  being  taken. 
An  approximately  cylindrical  volume  0.001-inch  in  diameter  and 
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0. 010-inch  in  length  was  chosen  as  best  meeting  the  previously  stated 
requirements,  allowing  statistically  significant  particle  counts  with 
the  chosen  particle  densities,  and  allowing  sharp  optical  definition. 
As  shown  in  Figure  4,  two  such  volumes  0.010-inch  apart,  aligned 
with  the  principal  flow  direction,  are  used  to  time  those  particles 
which  intercept  both  volumes.  The  average  transit  time  of  these 
during  the  period  of  measurement  is  taken  as  the  flow  velocity.  All 
particles  passing  through  one  of  the  volumes  are  counted  to  obtain 
the  particle  count.'  A  short  time  period,^  ,  on  the  order  of  0.01 
seconds  is  used  for  each  such  count  with  its  accompanying  velocity 
measurements.  A  long  time  period,  T,  on  the  order  of  30  seconds,  is 
used  for  the  complete  measurement  at  given  point  in  the  flow.  Each 
T  period  is  made  up  of  several  thousand  T  periods,  each  representing 
an  individual  concentration  measurement.  All  data  are  automatically 
recorded  as  they  are  taken  on  digital  tape  in  suitable  form  for 
direct  computer  processing. 
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The  concentration  calculated  from  a  single  %  period  is 


c  =  -EL 


where 


A  *  U 


U  = 


IU 


N 


The  mean  concentration  is  then 
7T  -  *C  , 

The  mean  velocity,  for  this  purpose, 
IU 

0  =  Nc 


(15) 

(16) 

(17) 


(18) 
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and  the  concentration  and  velocity  variations  for  a  given  T  period 
ore 


c  = 

C  -  C 

(IS) 

u  - 

# 

u  -  u 

(20) 

The  mean  product 

of  these  is 

w  Icu 

c*u  =  N<j- 

(21) 

The  product  which  is  needed,  however,  is  cv  •  Since  v 
cannot  be  measured  with  the  present  instrumentation,  it  is  necessary 
to  turn  to  the  literature  for  a  way  out  of  this  difficulty. 


Both  n^^and  have  been  measured  for  channel 
Therefore,  for  the  present,  since  c  is  a  scalar  quantity, 
proposed  that  the  following  relationship  be  used: 

cv  =  ciu 


Ttf 


f low(^). 
it  is 


(22) 


An  improved  version  of  the  instrument,  which  will  be  able  to 
measure  v  directly,  is  already  being  planned. 


The  instrument  as  constructed  consists  of:  (1)  traversing 
and  mounting  system  made  up  of  a  72-inch  lathe  bed  with  a  traveling 
carriage  permitting  the  sensor  system  40  inches  of  movement  along 
the  bed  and  6  inches  in  each  of  the  other  two  coordinate  directions; 
(2)  a  40-inch  long,  6-inch  square,  glass-walled  test  section  which 
mounts  over  the  lathe  bed;  (3)  sensor  system  arranged  physically 
as  shown  in  Figure  3  and  schematically  as  shown  in  Figure  5  (Two 
such  systems  arranged  to  move  independently  of  each  other  are 
planned  for  an  improved  version  of  the  instrument.);  (4)  electronics 
to  count  and  time  the  particles  and  put  the  information  in  proper 
'  format  for  recording;  (5)  digital  tape  recorder  to  store  the 
.  information. 
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Once  the  instrument  has  been  adjusted  to  scan  a  chosen  point 
in  the  flow  stream,  the  sequence  of  events  in  making  a  measurement  is 
as  follows:  (1)  The  operator  places  an  End-of-File  character  on  the 
tape  to  indicate  the  start  of  a  measurement. 

(2)  The  operator  places  a  six-character  label  on  the 
tape  to  identify  the  measurement.  This  step  may  be  repeated  as  many 
times  as  desired  to  put  additional  instructions  on  the  tape  con¬ 
cerning  the  particular  measurement.  This  information  must  be 
numerical . 

(3)  The  operator  starts  the  measurement  period  (T). 

(4j)  The  instrument  starts  the  first  counting  period(£j). 

(aj)  The  first  particle  velocity  is  recorded.  The 
following  information  is  given: 

J_  A  single-character  sensor  system  identifier. 

2  A  four-character  transit  time  measurement 
(50  nanosecond  resolution). 

2  A  nine- character  measurement  of  time  from 
the  start  of  the  counting  period  (2  microsecond  resolution). 

(a2)  The  second  particle  velocity  is  recorded. 


(an)  The  particle  velocity  is  recorded. 

(b)  The  counting  period  ends. 

(c)  A  six-character  counting  period  identifier 

is  recorded. 

(d)  A  six-character  total  particle  count  for  the 
counting  period  is  recorded. 

(e)  An  End-of-Information  character  is  recorded 
to  mark  the  end  of  the  recording  for  that  counting  period. 

•*  (42)  The- instrument  starts  the  second  counting 

period  (t^). 


(4n)  The  instrument  starts  the  Nth  counting  period  (TN), 
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(5)  The  measurement  period  ends.  The  sensor  system 
can  then  be  moved  to  scan  the  next  point  where  a  measurement  is  to 
be  made. 


When  a  series  of  measurements  have  been  completed,  the  taped 
information  is  fed  directly  into  a  digital  computer  programmed  to 
compute  the  desired  quantities. 

While  the  instrument  was  designed  to  meet  the  requirements 
of  the  problem  discussed  in  this  paper,  it  is  evident  it  has  wider 
applicability.  The  particles  in  the  flow,  themselves  of  concern  here, 
also  might  in  other  problems  serve  solely  as  tracers  for  measuring 
the  turbulent  characteristics  of  the  flow.  Since  no  physical  probe 
is  used,  investigation  of  flow  characteristics  is  possible  in  cases 
which  previously  were  difficult  or  impossible. 

The  improved  version  of  the  instrument,  which  is  now  being 
planned,  will  incorporate  two  independent  sensor  systems,  a  more 
accurate  and  flexible  traversing  system  to  carry  them,  particle  size 
discrimination,  the  ability  to  measure  transverse  velocity  components, 
and  such  improvements  in  electronics  as  are  available  at  the  time  it 
is  constructed.  This  will  make  possible  direct  measurement  of  the 
transverse  velocity  variation  -  concentration  variation  correlation, 
measurement  of  the  microscale  of  turbulence,  and  measurement  of  the 
correlation  coefficients  for  the  flow  itself. 

I 

It  is  hoped  that  significant  improvements  in  our  under* 
standing  of  turbulent  flow  can  be  made  through  the  use  of  this 
development  in  instrumentation. 


Figure  6 


206 


NOMENCLATURE 
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A  Cross  sectional  area  of  sensitive  volume  perpendicular 

to  main  flow  (equation  15) 

A  Constant  (equation  9,  11»  14) 

B  Constant  (equation  9,  11,  14) 

C  Constant  (equation  9,  11,  14) 

C  Concentration  (equations  1,  2,  3,  4,  15,  19) 

(T  =  C  (£,7),  Mean  concentration 

(T0  =  C  (0,0) 


JSt 
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N. 


Ntr 

R* 

Sc 

T 

U 

ui 


U*  Vr. 

V  * 


Molecular  diffusion  coefficient  for  particles 
Constant 

Number  of  particles  counted  during  period  T 
Number  of  particles  timed  during  period  T 
Number  of  periods  in  T 
Wal 1 -friction  velocity  Reynolds  Number 
Schmidt  number 

Time  for  complete  measurement  at  a  point 
Velocity 

Notation  representing  velocities  in  the  three 
coordinate  directions 

Average  velocity  of  particles  timed  during  T 

Velocity  of  an  individual  particle 

Wal 1 -friction  velocity 

Mean  velocity  in  the  x-direction 

Notation  representing  the  mean  velocities  in  the 
three  coordinate  directions 

Mean  velocity  in  the  y -direction 
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Mean  velocity  in  the  z-direction 
Dimensionless  mean  velocity  in  the  x-direction 
Variation  from  the  mean  concentration 


Notation  representing  the  mean  products  of  the 
concentration  variation  with  the  velocity  variations 
in  the  three  coordinate  directions 

Mean  product  of  the  concentration  variation  with 
the  x-direction  velocity  variation 

Mean  product  of  the  concentration  variation  with 
the  y-di recti on  velocity  variation 

Dimensionless  correlation  coefficient  betv/een  the 
concentration  variation  and  the  y-direction 
velocity  variation 

Integer  constant 

Time 

Variation  from  mean  velocity  in  the  x-direction 
2 

Mean  value  of  u 

Variation  from  mean  velocity  in  the  y-direction 
2 

Mean  value  of  v 

Coordinate  of  main  flow  direction 

Notation  representing  the  three  coordinate  directions 

Incremental  length  in  the  x-direction 

Coordinate  perpendicular  to  the  channel  walls 

Half  width  of  the  channel 

Incremental  length  in  the  y-direction 

Coordinate  perpendicular  to  the  main  flow  direction 
and  parallel  to  the  channel  walls 

Incremental  length  in  the  z-direction 

Dimensionless  y-direction  coordinate 


\  ' 
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*=  , 
M- 


S  * 


4  -  4(7) 
T=Y(t.  7) 


% 


Function  defined  by  equation  (13) 
Absolute  fluid  viscosity 
Kinematic  fluid  viscosity 
Dimensionless  x-direction  coordinate 
Fluid  density 
Shear  stress  at  the  wall 
Time  for  single  counting  period 
Function  defined  by  equetion  (14) 
Function  defined  by  equation  (10) 
Dimensionless  mean  concentration 

Particular  solution  to  equation  (6) 
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